A New Probe of the Molecular Gas in Galaxies: 
Application to M101 
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ABSTRACT 

Recent studies of nearby spiral galaxies suggest that photodissociation regions (PDRs) are 
capable of producing much of the observed HI in galaxy disks. In that case, measurements of the 
observed HI column density and the far-ultraviolet (FUV) photon flux responsible for the pho- 
todissociation process provide a new probe of the volume density of the local underlying molecular 
hydrogen. We develop the method and apply it to the giant Scd spiral M101 (NGC 5457). The 
HI column density and amount of FUV emission have been measured for a sample of 35 candi- 
date PDRs located throughout the disk of M101 using the Very Large Array and the Ultraviolet 
Imaging Telescope. We find that, after correction for the best-estimate gradient of metallicity in 
the ISM of M101 and for the extinction of the ultraviolet emission, molecular gas with a narrow 
range of density from 30-1000 cm -3 is found near star-forming regions at all radii in the disk of 
M101 out to a distance of 12' w 26 kpc, close to the photometric limit of R25 ~ 13.5'. 

In this picture, the ISM is virtually all molecular in the inner parts of M101. The strong 
decrease of the HI column density in the inner disk of the galaxy at Rq < 10 kpc is a consequence 
of a strong increase in the dust-to-gas ratio there, resulting in an increase of the H2 formation 
rate on grains and a corresponding disappearance of hydrogen in its atomic form. 

Subject headings: galaxies: individual (M101) — galaxies: ISM — ISM: clouds — ISM: molecules — HII 
regions — radio lines: galaxies — ultraviolet: galaxies 



1. Introduction 

Massive gas clouds in the ISM are considered to 
be the progenitors of young stars in the conven- 
tional picture of star formation in galaxy disks. 
The gas may initially be either in an atomic 
or in a molecular state (Elmegreen 1995) de- 
pending on the ambient physical conditions (UV 
flux and gas pressure), but as the volume den- 



sities increase through self-gravity, the gas be- 
comes mostly molecular on its way to forming 
stars. Once the gas is molecular and at densi- 
ties above « 100 cm~ 3 , the cooling rate is so high 
(e.g. Goldsmith & Langer (1978)) that rapid col- 
lapse follows. The rate-determining step in the 
conventional picture is then the aggregation of HI 
or diffuse H2 and the formation of Giant Molec- 
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ular Clouds. Elmegreen (1993) has discussed this 
problem in detail. 

For many late-type systems, the rate-determining 
step along the path to star formation appears to be 
the conversion of HI clouds into H 2 . The fraction 
of the ISM in atomic form is thought to increase 
as one progresses along the Hubble Sequence to- 
wards the late-type galaxies (see e.g. the review 
by Roberts & Haynes (1994)). Furthermore, the 
atomic fraction is thought to dominate the molec- 
ular by an order of magnitude or more for a sig- 
nificant number of actively-star-forming galaxies 
of types Sc and later (Young & Knezek 1989). 

However, in recent years evidence has been ac- 
cumulating that a significant fraction of the HI in 
the ISM of galaxies is in fact not a precursor to the 
star formation process, but instead is produced by 
photodissociation of the H 2 by UV photons em- 
anating from nearby newly-formed massive stars. 
In that case, either the way we compute the molec- 
ular fraction in galaxies via the CO(l-O) emission 
is wrong, or our ideas of how to control the other- 
wise runaway collapse of molecular clouds to form 
stars need further modification. Shu (1997) has 
summarized the view that magnetic fields may be 
a controlling factor in cloud collapse. 

The dissociation process was first described by 
Stecher & Williams (1967). The H 2 absorbs pho- 
tons emitted at 1108 A and 1008 A via electronic 
transitions in the Lyman (X 1 ^ — ► i? 1 E+) and 
Werner (X^j -» C 1 ^) bands. In the sub- 
sequent decay to the vibrationally excited lev- 
els of the ground electronic state, ~ 10% of the 
H 2 molecules will dissociate. Photons with wave- 
lengths as long as w 1850 A can continue to create 
HI by dissociating the "pumped" (X 1 ^, < v < 
14) H 2 via additional Lyman and Werner band 
transitions. The UV fluorescence spectrum pre- 
dicted from this process was first observed by Witt 
et al. (1989) in the nebula IC 63. 

The initial evidence that HI may be a product of 
the star formation process rather than a precursor 
to it was discovered in a comparison of the mor- 
phology of dust lanes and HI ridges in a prominent 
spiral arm of the galaxy M83 (NGC 5236) by Allen 
ct al. (1985, 1986). This comparison showed that 
the ridges defining the HI spiral arms peaked not 
on the dust lanes, as would have been expected in 
the conventional picture of density-wave-triggered 
star formation (Roberts 1969), but instead they 



were co- linear with the chain of HII regions created 
by young massive stars forming downstream from 
the dust lane. Other studies followed on M83 and 
other nearby spiral galaxies and have generally 
reached the same conclusions: for M83 (Tilanus & 
Allen 1993) , for M51 (Vogel et al. 1988; Tilanus 
& Allen 1987; Tilanus ct al. 1988; Tilanus & Allen 
1989, 1990, 1991; Rand et al. 1992; Knapen et al. 
1992), and for M100 (Knapen & Bcckman 1994, 
1996). 

Infrared data are providing important sup- 
port for this picture. For instance, models con- 
strained by KAO measurements of the 158/xm [CII] 
line indicate that photodissociation may produce 
as much as 70%-80% of the HI in the disk of 
NGC 6946 (Madden et al. 1993). More recently, 
spectra of nearby galaxies taken with the ISO 
satellite indicate that the bulk of the mid-infrared 
emission from galaxy disks arises in photodissocia- 
tion regions (PDRs) (Laurent et al. 1999; Roussel, 
Sauvage, & Vigroux 1999; Vigroux et al. 1999). 

Photodissociation therefore appears to be an 
important, and in some cases perhaps even the 
dominant, factor in defining the appearance and 
distribution of HI in galaxies, at least over the 
main body of their stellar disks. In the PDR pic- 
ture, HI is produced on the surfaces of molecu- 
lar clouds by the action of far-UV photons from 
nearby O and B (and some A) stars. The physics 
of PDRs has been worked out in detail (see e.g. the 
review by Hollenbach & Tielens (1999)), partly in 
response to new observational results on the in- 
frared spectra of nearby Galactic star-forming re- 
gions such as Orion. 

Recognition of the layered structure of the dif- 
ferent atomic and molecular species which exist 
in dense, bright PDRs like the Orion Bar (n « 

5 x 10 4 cm~ 3 , Tielens et al. (1993)) requires linear 
resolution better than ~ 0.02 pc, and the HI layer 
in this case would be only « 0.015 pc thick. How- 
ever, PDRs are also formed by modest UV fluxes 
incident on low-density clouds, and the HI layers in 
this case are much thicker. For example, Williams 

6 Maddalena (1996) have identified a large, low- 
density PDR in the Galaxy, near the molecular 
cloud G216-2.5. In this case, the HI layer is prob- 
ably 10-50 pc thick (n s=s 50 — 10 cm -3 ), and ap- 
pears in projection as an object about 50 x 200 
pc in size. While the IR emission from such large, 
low-density PDRs is hard to detect (but see Luh- 
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man & Jaffe (1996)), the HI emission is actually 
easier to measure, since the PDR is larger and 
therefore covers a larger fraction of the radio tele- 
scope beam at A21 cm. The problem with such 
measurements in the Galaxy which Williams and 
Maddelena had to circumvent is the confusion of 
the HI profiles by unrelated gas along the same 
line of sight. Synthesis imaging radio telescopes 
can now offer linear resolutions of order 100-200 
pc in the nearby galaxies, and the confusion prob- 
lem can be minimized by choosing galaxies which 
are viewed more face-on. 

The combination of high spatial resolution (~ 
100—200 pc) HI mapping along with far ultraviolet 
(FUV) and Ha imagery provides a unique oppor- 
tunity to explore the details of the photodissocia- 
tion picture in external galaxies. FUV imagery at 
w 1500 A (« 8.3 cV) directly measures the lower 
energy photons which can dissociate "pumped" 
molecular gas, and serves as a tracer for the more 
energetic photons in the range 912 < A < 1108 
A which start the pumping process for ground- 
state H 2 . Ha imagery is used to assess the impor- 
tance of extinction, which may be substantial in 
the far ultraviolet. Allen et al. (1997), hereafter 
A97, have used this technique with the Sb(r)I-II 
galaxy M81, finding that patches of HI emission 
are often located at the periphery of star form- 
ing regions mapped in the FUV, as expected for 
PDRs. The HI column densities in M81 are con- 
sistent with the photodissociation of a low density 
molecular gas by relatively modest UV fluxes 1 . 

In this paper, we extend the work of A97 to the 
Sc(s)I galaxy M101 and carry out a much more ex- 
tensive and quantitative analysis of the data. Ow- 
ing to its proximity (7.4 Mpc, Kelson (1996)), its 
nearly face-on orientation, and the availability of 
FUV, Ha, and HI data, M101 is the next obvious 
candidate for this type of work. The best angu- 
lar resolution presently available in HI for M101 is 

6"; this corresponds to w 215 pc, just adequate 
to permit identification of the larger PDRs, but 
still not sufficient to provide much morphological 
detail. The M101 data are discussed in §2; they 
have all been drawn from existing archival data 
sets. A representative sample of PDRs is identi- 

However, A97 did not adjust their results for the possible 
effects of internal extinction of the UV flux in M81, nor for 
changes in the dust-to-gas ratio. As we shall see in this 
paper, both corrections can be important. 



fied and discussed in §3. The volume density of 
the gas, n, associated with each PDR is derived 
in §3.5. The implications of the resulting radial 
profile of n are discussed in §4. 

The method used in this paper provides a new 
probe for studying the molecular gas content of 
galaxies which can be compared with the conven- 
tional results obtained from using the CO molec- 
ular emission (e.g. Solomon ct al. (1983)). Ques- 
tions have been raised about the reliability of CO 
as a tracer for molecular hydrogen in galaxy disks 
(e.g. Allen (1996); Combes (1999)), and it is of 
interest to explore other ways of establishing the 
amount and physical state of the ISM in different 
parts of galaxies. 

2. Archival Data Used 

The FUV image of M101 was obtained by the 
Ultraviolet Imaging Telescope (UIT) (Stecher et 
al. 1992, 1997) during the Astro-2 mission with 
the Bl filter (A = 1520 A, AA = 354 A) and an 
exposure time of 1310 s. The spatial resolution 
and pixel spacing of the image are 3" FWHM and 
1.13", respectively. Astrometric and photometric 
uncertainties are a — 1.8" and ~ 15%, respec- 
tively. An image of the Ha emission with resolu- 
tion « 2" FWHM was obtained from L. van Zee 
and M. Haynes. The VLA 21 cm line data were 
provided by R. Braun. The spatial resolution is 6" 
FWHM (215 pc), and the pixel spacing is 2" (72 
pc). The velocity channels are separated by 5.16 
km s _1 and the velocity resolution is 6.2 km s _1 . 
We re-derived a map of the HI column density, 
A (HI), by using the pixels in the channel maps 
whose values exceed 3 times the la noise value of 
16.8 K. The radial profile of the HI column den- 
sity derived from our data agrees with the original 
data portrayed in Figure 8a of Braun (1997). Full- 
size images of the FUV, Ha, and HI emission from 
M101 can be found in Waller et al. (1997), Van Zee 
ct al. (1998), and Braun (1995), respectively. 

The original astrometry of the FUV and Ha im- 
ages appears to be excellent based upon the gen- 
eral agreement in the positions of FUV peaks and 
Ha sources, and in the large-scale spatial distri- 
bution of the FUV continuum and the Ha and 
21 cm line emission. The astrometry of the HI 
images is intrinsically < 1", owing to the calibra- 
tion procedures used at the VLA. The FUV and 
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Ha images are then smoothed and resampled to 
match the spatial resolution and pixel spacing of 
the HI data. 

3. Results and Analysis 

The relationship between the FUV, Ha, and HI 
emission is illustrated in Figures 1 and 2. Boxes 
indicate the locations of regions shown in more 
detail in Figures 3 through 6. These regions are 
chosen arbitrarily to provide examples of the de- 
tailed FUV, Ha, and HI morphology. A large- 
scale spatial correspondence between the 21 cm 
line emission and the FUV and Ha emission is ob- 
served, as expected in both the conventional and 
photodissociation pictures. Star formation occurs 
where "sufficient amounts" of gas have accumu- 
lated in the conventional picture. However, that 
picture does not yet provide a quantitative con- 
nection between e.g. the amount of precursor gas 
(both atomic and molecular) which is present and 
the number of young stars formed from it. Deter- 
mining this relation from the observations is cur- 
rently an active area of research in the study of 
star formation in galaxies (e.g. Kennicutt (1998)). 

On the other hand, the photodissociation pic- 
ture provides a direct, quantitative connection be- 
tween the FUV flux emanating from a region of 
recent star formation, the volume density of the 
H 2 , and the column density of the surrounding 
HI. We have chosen in what follows to interpret 
our results in terms of the photodissociation pic- 
ture. While we clearly favor such an approach, 
we must emphasize that our data, by themselves, 
do not provide a definitive argument against the 
conventional picture. 

3.1. FUV and Ha: Mapping Star Forming 
Regions 

Extinction effects are much stronger at 1500 
A than in the visual, as indicated by reddening 
curves for the Milky Way (e.g. Ai 50 o = 2.6Ay 
(Savage & Mathis 1979)) and the Small Magellanic 
Cloud (e.g. A1500 = 5.6A V (Hutchings 1982)). 
This fact raises the concern that ultraviolet-based 
studies may be biased if entire star forming re- 
gions are totally obscured at FUV wavelengths. 
For this reason, A97 explore the effects of extinc- 
tion in M81 by comparing the FUV and Ha emis- 
sion. These authors find a strong correspondence 



between the small-scale structure of the FUV and 
Ha emission in M81, with every bright, reliable Ha 
source having a FUV counterpart. This suggests 
that the FUV morphology is relatively unaffected 
by extinction, or at least that the FUV and the 
Ha are affected in the same way. 

To explore the importance of extinction in 
M101, we examined a set of non-stellar Ha sources 
for FUV counterparts. We selected objects with 
the simple goal of constructing a set of represen- 
tative sources distributed across the galaxy which 
could be examined interactively for FUV coun- 
terparts. The automated point source extraction 
routines found in GIPSY and IRAF are sufficient 
for this purpose, yielding « 150 Ha sources when 
using a peak intensity threshold of 15o\ The re- 
sulting set does not form a statistically complete 
sample. Visual inspection of the data reveals that 
approximately 10% of the sources lack a detectable 
FUV counterpart within a radius of 10" (360 pc), 
at the present sensitivity. A comparison with 
the Ha images of M101 presented in Hodge et 
al. (1990) shows that these sources are in fact 
real HII regions, and not artifacts. Examples of 
Ha sources without FUV counterparts may be 
seen in Figure 3; these sources tend to be fainter 
Ha peaks. Source coordinates are given in Ta- 
ble 1. Sources such as those observed at RA = 
14' l 2 m 9.96 s , Dec= +54°35'16.5", and at RA = 
14^2™0.71 s , Dec = +54°32'2.7" arc foreground 
stars, based upon the continuum and narrow-band 
imagery of Hodge et al. (1990). The existence of 
equally faint Ha sources with FUV counterparts 
indicates that the lack of FUV emission is more 
likely to be an extinction effect than a sensitivity 
issue. 

We conclude that FUV wavelengths effectively 
map the locations of the majority of the star form- 
ing regions in M101, without severe extinction ef- 
fects. (The amount of extinction will be discussed 
further in Section 3.5.2). In this case, FUV map- 
ping offers an additional advantage over Ha map- 
ping in that FUV wavelengths will trace older sites 
of star formation in which ionizing stars are no 
longer present. Such FUV-bright/Ha-dim sites 
are powered primarily by B stars, whose photon 
energies are still favorable for the dissociation of 
H 2 . A97 find that 10% of 144 FUV sources in 
M81 do not exhibit Ha emission. We identified 
w 350 FUV sources distributed across M101, us- 
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ing a threshold peak intensity of 3a. Of these, 
the brightest 177 (half of the sources) were visu- 
ally examined for the presence of Ha counterparts. 
Approximately 7% of the 177 sources do not ap- 
pear to have an Ha counterpart within 10" (360 
pc), at the present sensitivity. Examples of such 
sources are seen in Figure 5 at A) RA(1950) = 
14 h m 40.28 s , Dec(1950) = +54°31'42.1"; B) RA 
= 14 h l m 0.62 s , Dec = +54°31'53.7"; C) RA = 
U h l m 3A5 s , Dec = +54°32'45.3". 

In the remainder of the paper, we focus our at- 
tention on a smaller sample of approximately 100 
FUV sources. The selection criteria were driven by 
our desire to identify a set of representative sources 
distributed across M101, manageable in size. The 
resulting sample is not intended to be statistically 
complete. The sample contains all FUV sources 
with peak intensities > 10tr identified by the point 
source extraction routines. These sources com- 
prise 2/3 of the set. The remaining sources were 
chosen by eye to bring the total number of sources 
to w 100, while keeping the distribution uniform 
across the galaxy. This technique ensures that we 
do not exclude sources that may have been missed 
by the point source extractors. The added sources 
have peak intensities of 3 — 10cr. All but one (a 
weak elongated source) belong to the set of 350 
sources discussed above. The properties of the 
Ha and HI emission are not considered during the 
source selection. 

3.2. The Detailed Morphology of the FUV 
and HI 

The small-scale relationship between the FUV 
and the HI emission is shown in Figures 4 and 
6. The general distribution of the HI gas fol- 
lows that of the diffuse FUV emission, with en- 
hancements in the 21 cm line emission seen near 
FUV peaks, as in M81. All but 2 of the w 100 
FUV sources exhibit 21 cm line emission within 
800 pc of the source. Specific regions of interest 
in M101 include the area surrounding SN 1951, 
located at RA(1950) = 14 h 02 m 07.5 s , Dec(1950) 
= +54°36'04", at the inner edge of the promi- 
nent eastern spiral arm. Here, FUV emission is 
coincident with and extends from the site of the 
NGC 5462 HII region (RA(1950) = 14 h 02 m 06.8 s , 
Dec(1950) = +54°36'20", our Figure 4, see also 
Figure 7 of Viallefond et al. (1981)). Emission 
from the 21 cm line lies along the periphery of the 



HII region but is absent in the vicinity of SN1951 
itself. Within the framework of the photodissoci- 
ation picture, an explanation for this morphology 
could be that this is a region of intense past and 
present star formation, as evidenced by the FUV, 
Ha and radio continuum emission. Generations of 
supcrnovae prior to SN1951 may have evacuated 
much of the surrounding molecular gas, resulting 
in the apparent hole in the gas distribution. Non- 
ionizing photons may have dissociated the skins of 
the remaining molecular gas clouds, creating the 
observed HI which neatly curves around the exist- 
ing FUV-emitting regions. A sketch of the spatial 
relationships between the FUV, Ha, and HI emis- 
sion in this type of situation may be found in Case 
IV, Figure 8 of A97. 

A case in which the star forming regions may lie 
above the mid-plane of the galaxy is seen at RA 
= 14 h l m 55.68 s , Dec = +54°33'23.7" in Figure 4. 
Here, the NGC 5461 complex shows enhanced HI 
emission slightly offset from the main FUV source, 
as well as spatially coincident HI. The NGC 5447 
complex, with peaks at RA = 14 h m 43.07 s Dcc= 
+54°30'35.5" and RA = 14 h m 41.69 s and Dec = 
+54°30'47.4", in Figure 6 is another example of 
this morphology, which occurs when HI appears 
spread over the FUV-emitting region but without 
much associated extinction. The HI in this case is 
behind the FUV source along the line of sight, as 
in Case I, Figure 8 of A97. 

Another interesting morphological behavior is 
exemplified by the small FUV-bright/Ha-dim 
source located in the HI depression at RA = 
14 h m 40.30 s Dec = +54°31'41.4", just north of 
NGC 5447 in Figure 6. This site may represent an 
older star forming region which has blown a hole 
through the plane of the galaxy. Non-ionizing 
photons from the remaining B stars could be dis- 
sociating the walls of the hole, producing the 
observed HI. This scenario corresponds to that 
illustrated in Case III, Figure 8 of A97. 

Finally, we note that the large loop of HI cen- 
tered at approximately RA = 14 /l l m 42 s , Dec = 
+54°31'45" in Figure 4 presents a case where some 
of the PDRs may be obscured. The northern por- 
tion of the loop shows the presence of fainter Ha 
sources and only weak FUV emission. The south- 
ern portion of the loop contains multiple, brighter 
FUV sources, however. Some of these sources also 
appear to lack Ha counterparts, and may be older 
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sites of star formation (e.g. the FUV peak located 
at RA = U h l m 39.12 s , Dec = +54°31'31.9"). This 
behavior could be explained if star forming regions 
in the northern portion of the loop lie closer to the 
mid-plane of the galaxy, in a behavior intermedi- 
ate between Cases I and II in Figure 8 of A97. 
The northern sources may also be younger sites of 
star formation which have not had sufficient time 
to break through the ISM. These types of sources 
are clearly in the minority, however; the locations 
of the majority of candidate PDRs are effectively 
mapped by the FUV emission. 

3.3. Candidate PDRs 

We have used the UIT data to select a group of 
candidate PDRs for further study. Of the « 100 
FUV-bright sources, we identified 35 regions for 
which 1) the total FUV flux and 2) the distance 
between the FUV source and the peaks of the 
associated 21 cm line emission could be reliably 
measured. The selected sources are therefore rel- 
atively isolated and in areas containing detectable 
HI. The primary selection criterion is in fact con- 
fusion, as 95% of the rejections were based upon 
the presence of another significant FUV source 
within < 15" (540 pc). Two sources were rejected 
due to a lack of 21 cm emission within a radius 
of « 15"; one additional source was rejected due 
to the weakness of the HI emission. The locations 
of the 35 selected regions are shown in Figure 7. 
The regions are classified into two subsets accord- 
ing to the intensity x °f the FUV flux seen by the 
neighboring HI clouds. The HI column density is 
directly related to x m the photodissociation pic- 
ture (see Section 3.4). Regions in Figure 7 marked 
with small square boxes are characterized by a nar- 
row range of \ such that 0.9 < \ < 10- These re- 
gions represent a set of "standard UV fluxes." Re- 
gions designated by open circles have FUV fluxes 
outside of this range. The FUV, Ha, and HI mor- 
phologies of typical candidate PDRs are illustrated 
in Figures 8 through 9. 

The quantity \ is given by \ — Ffuv(7-4: x 
10 6 /Phi) 2 /F , where F Q = 2.64 x 10~ 6 ergs cm" 2 
s _1 A -1 is the solar neighborhood value at 1500 
A (A97, and references therein), Fpuv is the total 
FUV flux observed by UIT in ergs cm _2 s _1 A _1 , 
Phi is the distance between the FUV source and 
the peak of the neighboring HI emission in parsecs, 
and 7.4 x 10 6 is the adopted distance to M101 



(in parsecs). The value of Fq is based upon the 
values of the local interstellar radiation field given 
by Draine (1978) and illustrated in Figure 1 of Van 
Dishoeck & Black (1988). The solar neighborhood 
value varies weakly over the wavelength range of 
the UIT Bl filter; the adopted value of F would 
change by less than 10% if we were to convolve the 
Draine (1978) curve with the UIT Bl filter profile. 

The measurements of pni and Nhi, the asso- 
ciated column density, are dependent upon the 
morphology of the source. As seen in Figure 9, 
multiple patches of 21 cm emission neighbor each 
source. Each patch serves as an equally valid 
probe of the surrounding ISM. We determine pni 
by calculating the average HI surface brightness 
in successive annuli for each PDR. The annuli 
are centered upon the FUV source. The quan- 
tity Phi is defined as the radius at which the sur- 
face brightness profile first reaches a local maxi- 
mum. To measure the associated column density, 
Nhi, we identify the brightest patch of 21 cm line 
emission located at a distance pni from the FUV 
source. The adopted value of Nhi corresponds to 
the brightest pixel within the patch. The distance 
between the brightest pixel and the FUV source is 
always within 0.5 pixels of pni- This measurement 
of Nhi is preferred to the average value of Nhi 
obtained from the surface brightness profiles since 
the annuli include regions lacking 21 cm line emis- 
sion. The average values therefore underestimate 
the column densities in the local ISM. 

The measurements of Fpuv, Phi, X, an( i Nhi 
for each region are listed in Table 2 and displayed 
in Figures 10 through 13 respectively. Filled cir- 
cles represent the regions with 0.9 < x < 10- The 
values of the sky background and any surround- 
ing diffuse FUV emission have been removed from 
the measurements of Fpuv- This is accomplished 
by measuring the average FUV surface brightness 
in successive annuli centered on the FUV source. 
The amount of diffuse + sky background per unit 
area is measured at the radius at which the FUV 
surface brightness levels off in the resulting radial 
profile. This value is subtracted from the total 
FUV flux contained within this radius. The FUV 
fluxes (Figure 10) vary by a factor of 100, but 
show no clear radial trend. With regard to the 
correspondence between the FUV and HI emis- 
sion, Figures 9 and 11 show many cases in which 
the small-scale HI emission is slightly offset from 
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or surrounds a FUV peak. Roughly half (16/35) 
of the FUV sources display an enhancement in 
the HI emission within 250 pc. For another 16 
sources, the separation between the FUV and HI 
peaks ranges from 250 pc to ~ 500 pc. While the 
photodissociation scenario suggests that the HI en- 
hancements are associated with the FUV emission, 
we remind the reader that other interpretations 
of the data may be valid, and that the observed 
morphology may simply reflect the general distri- 
bution of HI enhancements in M101. The HI en- 
hancements for the remaining 3 sources are fairly 
distant from the FUV sources (phi > 500 pc) and 
may not be associated with the FUV source. The 
values of pjjj do not display a significant radial 
trend (Figure 11). The resulting values of \ cover 
a large range (Figure 12), with no obvious trends 
with radius in the data. Values of Nhi repre- 
sent the peak column density associated with each 
source. Figure 13 indicates that the values of Nhi 
cover a wide range and tend to increase with in- 
creasing radius. This situation reflects the central 
deficit of 21 cm line emission in M101, a feature 
which is commonly observed in galaxies. 

Finally, the projected separation between the 
FUV and Ha sources (pHa) is given (Figure 14) 
in order to further quantify the effects of extinc- 
tion discussed in §3.1. The majority of the FUV 
sources (85%, 30/35) have an Ha counterpart lo- 
cated within 6" (215 pc; the FWHM of the HI 
image), as seen in Figure 14. Of the remaining 
sources, 4 have Ha emission which is slightly off- 
set (215 pc < pHa < 315 pc). The remaining 
FUV source lacks an obvious nearby Ha counter- 
part, with the nearest Ha source located over 600 
pc away. This source may be an older star forming 
region. The separation between the FUV and Ha 
emission is independent of the distance from the 
center of M101. 

3.4. The Quantitative Link Between the 
FUV and the HI 

The column density of HI produced in a PDR 
is fundamentally linked to the amount of far- 
ultraviolet (FUV) emission incident upon it and 
the local molecular gas volume density. The pro- 
duction of HI is calculated with the same physics 
as the production of the H 2 near-infrared emis- 
sion lines by fluorescence. Following e.g. Sternberg 



(1988), 



where 



JV(HI) = - In 
a 



DG 



(i) 



AT(HI) = the HI column density, 

a — the effective grain absorption 

cross section per H nucleus, 
D = the average unattenuated H 2 

photodissociation rate, 
R = the H 2 formation rate coeffi- 
cient, 

X — the FUV intensity relative to 
the local ISRF, and 

n = the total volume density of the 
gas. 

The quantity G is a dimensionless function of a, 
the absorption self-shielding function /, and the 
column density of molecular hydrogen Nn 2 - 

G = I afe-^^AN^ constant. 



-/ 

Jo 



G becomes constant for large values of iVn2 due to 
self-shielding (Sternberg 1988). 

Our analysis ignores the ionized component of 
the HI, which is produced by photons of higher 
energy (> 13.6 eV) than those we measure (~ 8.3 
cV). If such photons are present in the flux im- 
pinging on the molecular cloud, then some of the 
HI would disappear into HII and we would be un- 
derestimating A^(HI). However, in general both 
the PDR models and the highest-resolution ob- 
servations in the Galaxy (e.g. figures 3 and 30 in 
Hollenbach & Tielens (1999)) show that most of 
the HII appears closer to the exciting stars than 
does the HI, so any correction ought to be small. 

Equation 1 is strongly dependent upon the 
dust-to-gas ratio (5 = Av/Nh) and weakly de- 
pendent upon the gas temperature (T), since 

1.883 x 10- 21 (5/S a ) cm- 2 , 
3.0 x 10" 18/ 
(a/a ) 1/2 G , 



a 
R 
G 



V/^T^VCn cmV 1 , 



where So and cto refer to values in the solar neigh- 
borhood, and Vf{T) represents the efficiency of H 2 
formation. The product T 1 / 2 j/i?(T) is likely to be 
constant to within a factor of 2 (Hollenbach et al. 
1971). 
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Equation 1 also contains a dependence upon the 
level of obscuration in the immediate vicinity of 
the FUV source. While the variable \ represents 
the intrinsic FUV flux associated with the star- 
forming region, we generally observe an attenu- 
ated FUV flux. Assuming any extinction associ- 
ated with the star-forming region is in the form 
of an overlying screen of optical depth t(FUV), 
X(observed) = X e~ T{FUV) . Since the ISM in the 
immediate vicinity of the FUV source will have 
been disturbed by stellar winds and any prior su- 
pernovae, we do not attempt to link t(FUV), the 
local obscuration of the FUV source at 1500 A, 
to that implied by 5 = Av/Nh, the larger-scale 
dust-to-gas ratio. 

Assuming solar neighborhood values of <jq = 
1.883 x 10~ 21 cm 2 , D = 5.43 x 10~ n s"\ R = 
3 x 10~ 17 cm 3 s" 1 , and G ~ 5 x 10~ 5 (Stern- 
berg 1988), and neglecting the weak temperature 
dependence of R, equation 1 becomes: 



AT (HI) 



5 x 10 
IS/So) 



20 



In 



90 X ( 5 
n \S 



-1/2 



1 



(2) 



where X = x(observed)e T ( FC/v ). 

The behavior of A 7 " (HI) as a function of X IS 
displayed for 5/5 = 0.2, t(FUV) = 0andn = 30, 
n = 300, and n = 3000 in Figure 15. (These 
values of S/S and t(FUV) are appropriate for 
the outer regions of M101, as discussed in Section 
3.5.) Values of AT (HI) > 10 22 cm" 2 arc not likely 
to be observed as the atomic gas probably becomes 
optically thick at this point: 



AT (HI) = 1.82 x 10 



18 



T s t(v)cIv 



-> 1.82 x 10 18 T s tAv 
« 10 22 cm" 2 , 

for spin temperatures of T s w 100K and profile 
FWHMs of Ai; « 20 km s" 1 typical of M101 
(Braun 1997), and for optical depths of t w 2.5, 
corresponding to a ratio between the brightness 
and kinetic temperatures of Tb/Tk ~ 0.9. This 
value is appropriate for the highest-brightness re- 
gions of M101, as indicated in Figure 8a of Braun 
(1997). Figure 15 also shows the measurements for 
each of the 35 candidate PDRs. The data indicate 
that the properties of observed regions in M101 are 
consistent with photodissociation of an underlying 
molecular gas of moderate volume density. 



3.5. The Volume Density 

Equations 1 and 2 also suggest that our mea- 
surements of the HI column density and the FUV 
emission may be used to estimate the molecular 
gas density in M101. Solving for the volume den- 
sity of the gas yields: 



n = 90*(- 



-1/2 



e AT HI (5/5o)/5xl0 20 _ 1 



(3) 



for N m = AT (HI) and X = x(observed)e T (' FW ). 

This method provides a new probe of the molec- 
ular gas volume density with several advantages, 
including: 1) the method is independent of the 
CO/H 2 conversion factor, 2) the FUV emission 
can be directly measured from the UIT data, 3) 
the 21 cm line emission is insensitive to the exci- 
tation conditions in the ISM. The main disadvan- 
tages are: 1) the results are strongly dependent on 
the value of 5 = Av/Nh, which will vary through- 
out the galaxy, 2) the relationships contain a weak 
dependence on the gas temperature which is not 
yet well understood, 3) the amount of extinction 
affecting the FUV fluxes must be evaluated, and 
4) the observed values of N(HI) represent upper 
limits since the unknown geometry of the individ- 
ual PDRs prohibits inclination corrections. This 
method is also biased towards regions where signif- 
icant amounts of gas are present, i.e. regions which 
are massive enough to have formed some dozens of 
O-B stars in the last 10-100 x 10 6 years. 

3.5.1. The Dust-to-Gas Ratio 

Values of n derived from the data in Table 2, 
assuming a constant Milky Way dust-to-gas ra- 
tio 5 /So = 1 and no extinction (t(FUV) = 0) are 
also given in Table 2 as n raw . These values suggest 
that the volume density of the gas may decrease as 
a function of radius. However, the metallicity in 
M101 also changes strongly with radius (e.g. Ken- 
nicutt & Garnett (1996) and references therein). 
Since metal-rich regions also tend to be dustier, 
the observed radial trend in volume density most 
likely reflects a change in the dust-to-gas ratio. 

To investigate this possibility, a measure of the 
dust-to-gas ratio is needed. Following Issa et al. 
(1990) and Schmidt & Boiler (1993), we assume 
that S = Av/Nh is directly proportional to the 
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metallicity, as traced by the oxygen abundance. 
Radial profiles of the quantity [12 + (logO/H)] for 
M101 are given in Kennicutt & Garnett (1996) 
for different values of the calibration between the 
R23 = ([Oil] + [O IIl])/H/3 abundance parameter 
and oxygen abundance. The radial profiles ap- 
pear to be linear for the Dopita & Evans (1986) 
and McCall et al. (1985) calibrations or to steepen 
somewhat in the central regions for the Edmunds 
& Pagel (1984) calibration. In the rest of the pa- 
per, we refer to the metallicity gradients in terms 
of the calibration source, but remind the reader 
that the gradients all refer to the Kennicutt & 
Garnett (1996) data. Figure 16 shows the range 
of radial dependence of the dust-to-gas ratio for 
M101. Values of n derived using dust-to-gas ra- 
tios based upon the Edmunds & Pagel type metal- 
licity gradient are shown in Figure 17. This figure 
shows values increasing from n ~ 10 cm~ 3 in the 
central star forming regions of M101 to n ~ 1000 
cm~ 3 in the periphery of the galaxy. The trend is 
similar for the Dopita & Evans or the MRS type 
metallicity gradients, although slightly lower val- 
ues are obtained (n ~ 0.1 cm~ 3 ton~ 100 cm -3 ). 
Values of n derived using an average of the Dopita 
& Evans and MRS type metallicity gradients are 
shown in Figure 18. 

3.5.2. Obscuration Effects 

The discussion up to this point has assumed 
that extinction effects in these FUV-bright re- 
gions are negligible, based upon the fact that the 
clear majority of Ha sources have FUV counter- 
parts. Studies of FUV-bright regions in M81 and 
M51 suggest that such regions are still obscured by 
dust with optical depths ranging from ry <~ 0.4 to 
t v ~ 3.2 (Hill ct al. 1995, 1997), however. These 
values would translate to non negligible redden- 
ings of at least 1 to 9 magnitudes at 1500 A. We 
use the Ha and H/3 measurements of Scowen et 
al. (1992) to apply an extinction correction to our 
values of \. These authors find that the central re- 
gions of M101 are characterized by higher optical 
depths than the outer regions, consistent with the 
general behavior of the dust-to-gas ratio. Based 
upon the scatter plot of optical depths in Figure 
3 of Scowen et al., we adopt a radial profile of the 
optical depth of the form ry = 2.3 — R/6.5, where 
R is the radius, for R < 15 kpc, and ry = for 
R > 15 kpc. These values provide a lower limit to 



the required extinction corrections, as radio obser- 
vations of HII regions in M101 indicate that opti- 
cal depths derived from optical Ha/H/3 line ratios 
only trace a fraction of the total optical depth in 
those regions (Viallcfond ct al. 1982). The en- 
suing extinction-corrected values of n are shown 
in Figure 19 for the Savage & Mathis (1979) ex- 
tinction curve (t(FUV) = 2.6ry) and Edmunds 
& Pagel (1984) type metallicity gradient. Even 
higher volume densities (approaching n ~ 2000 
cm~ 3 in the central regions) are obtained with 
the Kinney et al. (1994) extinction curve. Fig- 
ure 20 shows extinction-corrected values of n for 
the Savage & Mathis (1979) extinction curve and 
the average of the Dopita & Evans and MRS type 
metallicity gradients. 

These calculations use corrections based upon 
averaged extinction and dust-to-gas values; the 
volume densities could be further constrained by 
obtaining data on the specific optical depths and 
dust-to-gas ratios of the specific 35 candidate 
PDRs in our study. 

4. Discussion and Conclusions 

Figures 19 and 20 are our current best estimates 
for the radial dependence of the total gas volume 
density in M101. This gas should be essentially all 
H 2 . In the process of correcting the "raw" values 
of volume density (Table 2) for the radial varia- 
tions of dust-to-gas fraction (resulting in Figures 
17 and 18) and UV extinction (finally producing 
Figures 19 and 20), the radial variations have dis- 
appeared. The molecular gas shows a surprisingly 
narrow range of values from 30 - 1000 H2 molecules 
cm~ 3 for the Edmunds & Pagel type metallicity 
gradient, with no clear trend from the inner HI 
-deficient regions near Rq < 5 kpc, through the 
HI -rich "main body" of the galaxy near Rq ~ 15 
kpc, all the way out to the HI -poor outer parts at 
Rg > 25 kpc. We note that the Edmunds & Pagel 
type metallicity gradient is a good representation 
for the mean of an ensemble of nearby galaxies 
(Panagia 2000), which leads us to favor Figure 19 
as the most likely final result. 

One concern is that the narrow range of values 
in Figure 19 could result from observational selec- 
tion. For example, regions of high volume density 
located near typical FUV sources may be missed 
since the HI would have a relatively low column 
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density and a low filling factor. This is not an issue 
for this study since only 3% of our initial sample 
of FUV sources were rejected due to a lack of 21 
cm line emission. To confirm this statement, we 
examined an additional 10 distinct FUV sources. 
Measurable HI is found within 800 pc for all ten 
sources. This is consistent with the spatial corre- 
lation between FUV sources and 21 cm line emis- 
sion observed in Figures 2, 4, and 6. Virtually all 
FUV sources are associated with 21 cm line emis- 
sion, with the exception of the innermost regions 
of M101, where the HI column density is supressed 
due to the strong increase in the dust-to-gas ratio. 

What kind of observational selection could op- 
erate on the low-density side of Figure 19? Ac- 
cording to Equation 1, H2 with low volume den- 
sity can produce large columns of HI which may 
become optically thick, leading us by Equation 3 
to overestimate the H2 density and thus depop- 
ulate the lower part of Figure 19. However, gas 
of density e.g. n < 10 cm~ 3 with HI columns in 
the optically-thick regime of N(HI) w 10 22 cm~ 2 
would be typically > 300 pc in size. Examination 
of Figure 2 and the zoomed Figures 4 and 6 shows 
that the highest-brightness HI regions are unre- 
solved, however. NGC 5447 in Figure 6 is a good 
example of this, as are NGC 5461 and NGC 5462 
in Figure 4. This means the path lengths are 
< 200 pc, which implies values of N(HI) < 6 x 10 21 
cm~ 2 in Figure 15, where optical depth in the HI is 
not likely to be a problem. Therefore, the paucity 
of points at low volume densities in Figure 19 is 
very likely real. 

We conclude that the 35 regions studied do in- 
deed comprise a representative sample of PDRs, 
and that observational selection is not likely to be 
the cause of the narrow range of values for n(H 2 ) 
in Figure 19. 

Our results refer only to the ISM in the im- 
mediate vicinities of young stars and clusters of 
young stars. We do not know from the present 
observations if the results are representative of the 
ISM as a whole. Furthermore, the measurement 
of the HI emission does not immediately provide 
a total gas mass without additional assumptions 
about the geometry of the ISM since the HI emis- 
sion is a surface phenomenon in the PDR picture. 
Further discussion of the implications of our re- 
sults for the ISM content of galaxies in general is 
beyond the scope of this paper; we hope to return 



to these questions in a future publication. 

Waltcrbos & Braun (1996) have proposed an 
interpretation for the radial distribution of the 
bright HI peaks in M101 and other nearby spi- 
rals in terms of variations in the hydrostatic pres- 
sure of the ISM. Their picture is insightful and 
deserves a more thorough discussion in the light 
of the present results than space here permits. 

UV photons clearly dominate the physics of 
the ISM in the immediate neighborhoods of young 
stars. The production of HI from H2 by photodis- 
sociation is a natural and inevitable consequence 
of this physics. We have described how photodis- 
sociation can explain the geometrical structure of 
the HI on scales of ~ 100 pc in nearby galaxies 
and have discussed how the physics of photodis- 
sociation can be used to provide a new probe of 
the density of the underlying molecular hydrogen. 
In addition, we have presented a straightforward 
explanation for the disappearance of HI in the in- 
ner parts of galaxies. Acceptance of this picture 
requires a shift in parts of the paradigm for star 
formation from the ISM, from viewing HI as a pre- 
cursor, to seeing it as a product of the star forma- 
tion process. While the present results do not by 
themselves demand discarding the old paradigm, 
they further strengthen the viability of the new 
one. 
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paper. We are grateful to our referee for a care- 
ful reading of our manuscript and comments which 
improved the paper. The authors have made use of 
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which is operated by the Jet Propulsion Labora- 
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Fig. 1. — FUV and Ha emission. Contours of the FUV emission are superposed on a greyscale representation 
of the Ha emission. The contour levels are 1.5 x 1CT 17 , 5 x 1CT 17 , 1 x lCT 16 , 5 x 1(T 16 , and 1 x 10 -15 
ergs cm _2 s _1 A _1 . Boxes 1 and 2 indicate the regions shown in more detail in Figures 3 and 5. 



Fig. 2. — FUV and 21 cm line emission. Contours of the FUV emission are superposed on a greyscale 
representation of the 21 cm line emission. Contour levels and boxes are as in Figure 1. 



Fig. 3. — FUV and Ha emission - region 1. The detailed behavior of the Ha emission from region 1 (see 
Figure 1) is shown in grayscale. Contours of the FUV emission are given for FUV fluxes of 1.5 x 10 -15 , 
4 x 1(T 17 , 6 x 1(T 17 , 8 x 1(T 17 , 1 x 1(T 16 , 1.5 x 10" 16 , 2 x 1(T 16 , 3 x 1(T 16 , and 5 x 10~ 16 ergs cm-V 1 !" 1 . 
Ha sources without FUV emission are labelled according to their Hodge numbers. 



Fig. 4. — FUV continuum and 21 cm line emission - region 1. The detailed behavior of the 21 cm line 
emission from region 1 (see Figure 2) is shown in grayscale. Contours of the FUV emission are given as in 
Figure 3. The FUV contours appear smooth in comparison to the 21 cm line emission since the ultraviolet 
data were not resampled after smoothing; the FUV emission is essentially oversampled in this representation. 



Fig. 5. — FUV and Ha emission - region 2. The detailed behavior of the Ha emission from region 2 (see 
Figure 1) is shown in grayscale. Contours of the FUV emission are given for FUV fluxes of 1.5 x 1CU 15 , 
4 x 1(T 17 , 6 x 1(T 17 , 8 x 1CT 17 , 1 x 1CT 16 , 1.5 x 1CT 16 , 2 x 1CT 16 , and 5 x 1(T 16 ergs cm-V^" 1 . Sources 
A, B, and C are FUV sources lacking significant Ha emission. 



Fig. 6. — FUV continuum and 21 cm line emission - region 2. The detailed behavior of the 21 cm line 
emission from region 2 (see Figure 2) is shown in grayscale. Contours of the FUV emission are given as in 
Figure 5. 



Fig. 7. — FUV continuum emission. The FUV image shown here has been smoothed and regridded to match 
the spatial resolution and pixel spacing of the HI data. The locations of 35 relatively isolated sources are 
shown. Squares correspond to regions with 0.9 < \ < 10- Regions with \ values outside of this range are 
marked by circles. 
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Fig. 8. — PDR Morphology: FUV and Ha. This figure shows the FUV and Ha emission from eight of the 
candidate PDRs with 0.9 < x < 10. The regions arc in order of increasing distance from the nucleus of 
M101. The individual values of X are a) 5.8, b) 3.0, c) 3.9, d) 2.0, e) 2.3, f) 2.5, g) 7.0, and h) 1.8. The FUV 
contour levels are, in units of 1 x 10~ 16 ergs cm _2 s _1 A _1 , a) 0.25, 0.5, 0.75, 1, 1.5, 2, 2.5, and 3, b) 0.25, 
0.5, 0.75, 1, 1.5, 2, 2.5, 3, and 3.5, c) 0.4, 0.6, 1, 1.4, 1.8, 2.2, and 2.6, d) 0.25, 0.5, 1, 1.5, 2, 2.5, and 3, e) 
0.2, 0.4, 0.8, 1.2, 1.6, and 2, f) 0.2, 0.4, 0.8, 1.2, and 1.6, g) 0.25, 0.5, 1, 2, 3, 4, and 5, and h) 0.10, 0.20, 
0.30, and 0.40. 
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Fig. 8.— PDR Morphology: FUV and Ha, continued. 
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Fig. 9.— PDR Morphology: FUV contours and HI grayscale. As in Figure 8, for the FUV continuum and 
21 cm line emission. The HI is offset from and partially surrounds the FUV source. The area filling factor 
appears to be roughly ~ 30%, independent of radial distance. 
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Fig. 9. — PDR Morphology: FUV contours and HI grayscale, continued. 
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Fig. 10. — Observed FUV fluxes (Ffuv)- This figure illustrates the range of FUV fluxes produced by the 
35 isolated sources, as observed by UIT. Diffuse FUV emission, as well as any sky emission, have been 
subtracted. Filled circles correspond to PDRs with 0.9 < x < 10. PDRs with \ values outside of this range 
are denoted by open circles. The observed FUV fluxes do not show a radial trend. 
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Fig. 11. — Observed FUV-HI separation (pm)- The distance between each FUV peak and the peak of the 
surrounding 21 cm emission is shown. The quantity pm does not show a clear trend with radius. 
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Fig. 12. — Derived x Values. The FUV flux observed at the location of the peak HI is derived from Ffuv 
and phi- The values of x are independent of radius and clustered between x = 0.9 and x = 10. 
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Fig. 13. — Peak Column Densities (Nhi)- The peak column density in the vicinity of each FUV source is 
measured from the HI map. Column densities generally increase as the distance from the nucleus increases, 
reflecting the large-scale trend seen in Figure 2. 
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Fig. 14. — Observed FUV-Ha separation (pHa)- The distance between each FUV peak and the peak of the 
associated Ha emission is shown. 
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Fig. 15. — The Relationship between JV(HI) and \- The observed and predicted values of 7V(HI) are shown as 
a function of \- The modeled behavior of TV (HI) assumes values of 5/8q = 0.2 and t(FUV) — 0, appropriate 
for the outer regions of M101. The observations are clearly consistent with the physics underlying the 
photodissociation picture. 
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Fig. 16. — Dust-to-Gas Profiles. The illustrated dust-to-gas profiles are obtained from the metallicity 
gradients given in Kcnnicutt & Garnett (1996). The solid line is based upon the calibration of Edmunds & 
Pagel (1984); the dashed line reflects the average of the calibrations of Dopita & Evans (1986) and MRS. 
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Fig. 17. — Gas Volume Density with Varying Ay/Nn, but no Extinction Correction. In the case where the 
metallicity gradient steepens in the central regions (Edmunds & Pagel (1984) calibration), values of n range 
from n ~ 10 cm' 3 to n ~ 1000 cm~ 3 . The volume density appears to increase with radius, x has not been 
corrected for extinction. 
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Fig. 18. — As in Figure 17, but with an alternative metallicity correction using the average of the calibrations 
of Dopita & Evans (1986) and MRS. Note that % has not been corrected for extinction. 
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Fig. 19. — Gas Volume Density with varying Av/Nh, including an extinction correction. The volume 
densities have been derived using the Edmunds & Pagel type metallicity gradient and the Savage & Mathis 
reddening curve. (This is our preferred set of corrections at the present time.) 
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Fig. 20. — As in Figure 19, but with alternative mctallicity correction, and including an extinction correction. 
The volume densities have been derived using the average of the Dopita & Evans and MRS calibrations and 
the Savage & Mathis reddening curve. 
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Table 1 

A Sample of Ha sources lacking FUV counterparts 



Hodge Number 


RA (1950) 


Dec (1950) 


888 


14 h l m 38.5 s 


+54°38'34.4" 


891 


14 1 38.3 


+54 35 46.5 


939 


14 1 40.5 


+54 35 10.8 


979 


14 1 42.4 


+54 36 4.8 


1117 


14 1 56.9 


+54 33 44.2 


1122 


14 1 58.4 


+54 34 0.1 


1126 


14 1 58.9 


+54 34 7.0 



Table 2 

Observed and Derived Properties of the 35 Candidate PDRs. 



RA 


Doc 


Radius 


Ffuv 


Phi 


X 


Nhi 




PHa 




n b 


(1950) 


(1950) 


(kpc) 


(ergs cm~ 2 s~ 1 A _1 ) 


(pc) 




(atoms cm ) 


(PC) 


("m"™ 3 ) 


(cm" 3 ) 


14 h l m 18.77 s 


+54°35'26.9" 


2.5 


3.00 x 10~ 15 


359. 


0.48 


0.23 X 


10 21 


225. 


76. 


2047. 


14 1 13.93 


+54 36 24.9 


4.7 


1.21 x 10~ 15 


682. 


0.05 


0.53 x 


10 21 


110. 


2.6 


49. 


14 1 42.47 


+54 34 52.9 


5.0 


5.99 x 10~ 15 


144. 


6.03 


1.73 x 


10 21 


109. 


18. 


116. 


14 1 12.11 


+54 33 58.9 


5.3 


1.52 x 10~ 15 


144. 


1.53 


0.89 x 


10 21 


107. 


28. 


443. 


14 1 31.43 


+54 37 52.9 


5.9 


6.32 x 10~ 15 


430. 


0.71 


1.58 x 


10 21 


106. 


2.8 


33. 


14 1 26.37 


+54 32 16.9 


6.4 


14.27 x 10~ 15 


502. 


1.17 


0.97 x 


10 21 


168. 


18. 


322. 


14 1 08.87 


+54 36 50.9 


6.5 


7.21 x 10~ 15 


161. 


5.74 


1.45 x 


10 21 


147. 


30. 


482. 


14 1 47.07 


+54 34 32.8 


6.6 


9.45 x 10~ 15 


201. 


4.85 


1.73 x 


10 21 


147. 


14. 


216. 


14 1 13.24 


+54 37 48.9 


6.9 


4.58 x 10~ 15 


201. 


2.35 


1.15 x 


10 21 


167. 


24. 


454. 


14 1 10.04 


+54 32 50.9 


7.3 


18.47 x 10~ 15 


359. 


2.98 


2.61 x 


10 21 


314. 


1.5 


29. 


14 1 39.00 


+54 31 32.9 


8.9 


8.69 x 10~ 15 


215. 


3.89 


3.12 x 


10 21 


240. 


0.68 


39. 


14 58.07 


+54 34 30.7 


9.0 


3.32 x 10~ 15 


201. 


1.71 


1.15 x 


10 21 


110. 


17. 


429. 


14 1 55.57 


+54 33 24.7 


9.9 


31.81 x 10~ 15 


287. 


8.01 


5.86 x 


10 21 


118. 


0.0059 


4. 


14 1 37.89 


+54 39 50.9 


10.5 


15.98 x 10~ 15 


72. 


64.38 


3.00 x 


10 21 


120. 


15. 


1604. 


14 2 02.25 


+54 34 18.6 


11.3 


3.20 x 10~ 15 


144. 


3.22 


2.95 x 


10 21 


74. 


0.79 


106. 


14 1 56.54 


+54 38 10.7 


11.3 


3.48 x 10~ 15 


144. 


3.50 


3.70 x 


10 21 


126. 


0.19 


53. 


14 50.46 


+54 36 10.6 


11.4 


15.13 x 10~ 15 


395. 


2.02 


3.44 x 


10 21 


114. 


0.19 


42. 


14 1 55.40 


+54 38 44.7 


11.7 


1.46 x 10~ 15 


305. 


0.33 


2.29 x 


10 21 


144. 


0.30 


23. 


14 2 11.26 


+54 36 48.4 


14.3 


4.07 x 10~ 15 


108. 


7.28 


3.92 x 


10 21 


77. 


0.26 


120. 


14 1 14.89 


+54 28 46.9 


14.4 


38.71 x 10~ 15 


359. 


6.24 


3.48 x 


10 21 


106. 


0.53 


140. 


14 33.69 


+54 34 24.2 


16.6 


10.44 x 10~ 15 


305. 


2.33 


4.02 x 


10 21 


80. 


0.068 


53. 


14 34.20 


+54 32 12.2 


17.6 


16.23 x 10~ 15 


215. 


7.27 


4.23 x 


10 21 


90. 


0.14 


191. 


14 34.08 


+54 37 38.2 


17.2 


2.98 x 10~ 15 


108. 


5.35 


4.19 x 


10 21 


131. 


0.11 


127. 


14 25.65 


+54 34 02.0 


19.2 


2.34 x 10~ 15 


610. 


0.13 


2.29 x 


10 21 


139. 


0.12 


16. 


14 24.30 


+54 32 53.9 


20.1 


3.00 x 10~ 15 


72. 


12.10 


2.91 x 


10 21 


145. 


3.2 


1172. 


14 53.12 


+54 43 22.7 


20.4 


10.09 x 10~ 15 


287. 


2.54 


3.79 x 


10 21 


142. 


0.12 


164. 


14 2 25.61 


+54 39 38.1 


20.7 


27.62 x 10~ 15 


287. 


6.96 


3.40 x 


10 21 


111. 


0.70 


581. 


14 40.94 


+54 28 02.4 


21.1 


2.68 x 10 -15 


161. 


2.14 


2.69 x 


10 21 


100. 


0.88 


284. 


14 29.87 


+54 40 34.1 


21.1 


10.78 x 10" 15 


610. 


0.60 


1.11 X 


10 21 


640. 


6.7 


266. 


14 1 56.62 


+54 44 14.7 


21.5 


1.12 x 10~ 15 


466. 


0.11 


1.86 x 


10 21 


139. 


0.24 


26. 


14 22.74 


+54 30 49.9 


22.1 


1.51 x 10 -15 


395. 


0.20 


1.24 x 


10 21 


142. 


1.7 


91. 


14 2 30.28 


+54 41 57.9 


24.5 


2.30 x 10 -15 


161. 


1.83 


3.00 x 


10 21 


196. 


0.41 


392. 


14 1 50.63 


+54 46 06.8 


24.6 


3.93 x 10~ 15 


287. 


0.99 


3.00 x 


10 21 


188. 


0.22 


212. 


14 2 43.54 


+54 38 05.4 


24.7 


98.54 x 10~ 15 


251. 


32.41 


5.77 x 


10 21 


17. 


0.028 


2466. 


14 1 05.56 


+54 46 40.8 


25.5 


2.04 x 10~ 15 


251. 


0.67 


4.71 x 


10 21 


268. 


0.0049 


85. 



^Volume densities derived using the (unrealistic) case of Av/Nh — constant. 

b Volume densities derived using the Edmunds & Pagel metallicity gradient and Savage & Mathis extinction law. 
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